Coreless vortex dipoles and bubbles in phase-separated binary condensates 
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Vortex dipoles are generated when an obstacle moves through a superfluid. In case of phase- 
separated binary condensates, with appropriate interaction parameters in pan-cake shaped traps, 
we show that coreless vortex dipoles are created when a Gaussian obstacle beam traverses across 
them above a critical speed. As the obstacle passes through the inner component, it carries along a 
bubble of the outer component. Using Thomas-Fermi approximation, we show that phase-separated 
binary condensates can either support vortices with empty or filled cores. For time dependent 
obstacle potentials, ramped down in the present case, relative energy scales of the system influence 
the dynamical evolution of the binary condensate. 
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Vortices are among the most fundamental topological 
structures in fluid flows. These occur in a variety of flu- 
ids at different length scales and in different settings. In 
Bose- Einstein condensates (BECs), vortices carry inte- 
gral angular momenta and serve as evidence of super- 
fluidity. Formation and dynamics of solitary and multi- 
ple vortices in single species condensates is well studied. 
These vortices with empty cores and a phase singularity 
at the center of the cores are normal vortices, whose ana- 
logues are found in variety of other fluids. An important 
development in vortex dynamics is the realization of vor- 
tex dipoles, a pair of vortex and anti-vortex, due to super- 
fluid flow past an obstacle pQ. Furthermore, the recent 
development to observe vortex dipoles in real time [5] 
opens up the possibility to examine dynamics. The vor- 
tex dipoles can also be considered as the embedded sec- 
tion of a vortex ring, which were first observed in BEC as 
decay product of dark solitons [5J . Another phenomenon 
in superfluids, where vortex dipoles are crucial is tran- 
sition to quantum turbulence. The initial requirement 
is the formation of multiple vortex dipoles, followed by 
reconnections to form a vortex tangle, which then leads 
to quantum turbulence [4j. In multi-component conden- 
sates, coreless vortices [5] and Skyrmions [6 have been 
experimentally observed in spinor condensates. Core- 
less vortices in multi-component spinless condensates are, 
however, yet to be experimentally realized. 

Theoretically, the flow of the miscible binary and 
spinor condensates across a Gaussian obstacle potential 
has been investigated [7HS] ■ A related work is the seeding 
of vortex rings when a condensate bubble moves through 
the bulk of another condensate in binary condensates. 
This is examined in a recent work 10J. In addition, 
phase-separated binary condensates, under suitable con- 
ditions, are appropriate systems to observe phenomena 
related to dynamical instabilities. These include, modu- 
lational instability [11] [12] , Kelvin-Helmholtz instability 
[15] . Rayleigh- Taylor instability [21 [15], countersuper- 
flow instability [THl [T7] and Richtmyer-Meshkov instabil- 



ity [18]. Among these, the countersuperflow instability 
was recently realized experimentally [19] . Phase separa- 
tion, which has been observed experimentally [201 I2~I] - 
is also a necessary condition to observe coreless vortex 
dipoles with a moving obstacle. 

In this letter we report the theoretical study of vor- 
tex dipole formation when a Gaussian obstacle moves 
through phase-separated binary condensates. As a spe- 
cific case, we consider the binary condensate consisting 
of 85 Rb and 87 Rb. In this, one of the scattering lengths 
can be tuned by using magnetic Feshbach resonances [5D] . 
This allows us to choose values of scattering lengths suit- 
able for generation of coreless vortex dipoles. 



BINARY CONDENSATES AND VORTEX 
DIPOLES 

Dynamics of weakly interacting binary condensate at 
zero temperature is well described by a set of coupled GP 
equations 
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in mean field approximation, where i = 1, 2 is the species 
index. Here Un = Aitfi^an/mi, where rrii is the mass 
and an is the s-wave scattering length, is the intra- 
species interaction, Uij — 27rfi 2 a.y /my, where my = 
rriimj / (rrii+mj) is the reduced mass and is the inter- 
species scattering length, is the inter-species interaction, 
and Vi(r) is the trapping potential experienced by ith 
species. In the present work, we consider binary conden- 
sate consisting of 85 Rb and 87 Rb for which mi w rri2- 

Furthermore, we also consider identical trapping po- 
tentials, which are axially symmetric, for both the 
species. The total potential is then 
™, ,2 



V(T,t) 



(x 2 + y 2 + f3 2 z 2 ) + V obB (x,y,t), 



2 



where V obs (x,y,t) = V (t) exp{-2([x - x (t)} 2 + y 2 )/w 2 } 
is the blue detuned Gaussian obstacle potential and 
ft is the anisotropy parameter. Define the oscilla- 
tor length of the trapping potential a osc = y/H/ (rruJ), 
and consider hu as the unit of energy. We can 
then rewrite the equations in dimensionless form with 
transformations f = r/a osc , t = tu and (j)i(r,t) — 
yJa^ sc /Ni^!i(r, t). In pancake-shaped traps (/3 3> 1), 
0(r,f) = ip(x,y,t)C(z)exp(-il3t/2) [22], where C = 
(/3/(27r)) 1 / 4 exp(— /3z 2 /4) is the ground state wave func- 
tion in axial direction. The Eq. ([I]) can then be reduced 
to the two dimensional form 
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where ug = 2a it Ni V / 2Trj3. l /a osc and u tl = 
2ciij Nj y/2ir fti / a osc . Here we have neglected a con- 
stant term corresponding to energy along axial direction 
as it only shifts the energies and chemical potentials 
by a constant number without affecting the dynamics. 
In the present work, we consider U12 > ^/mh«22 so 
that the ground state of the binary condensate is 
phase-separated. Geometry of the density distribution 
is such that the species with the lower repulsion energy 
forms a core and the other species forms a shell around 
it. For convenience, we identify the former and later 
as the first and second species, respectively. With this 
labelling, interaction energies u\\ < U22 and for equal 
populations, this implies an < 022- 

To be specific, we consider 85 Rb- 87 Rb binary conden- 
sate with an = 460ao, 022 = 99ao, and ai2 = 214ao 
as the scattering length values and 2N\ — N2 = 10 6 as 
the number of atoms. Here an is tunable with magnetic 
Feshbach resonance [23]. With these set of parameters, 
the stationary state of 85 Rb- 87 Rb binary condensate is 
just phase-separated. The trapping potential and obsta- 
cle laser potential parameters are same as those consid- 
ered in Ref. [T], i.e. u/(2w) = 8Hz, a = 1, /3 = 11.25, 
Vb(0) = 93.0/kj, and wq — 10^m. Hereafter we term this 
set of scattering lengths, number of atoms and trapping 
potential parameters as set a. 

In hydrodynamics, the velocity field of a vortex dipole 
is the vector sum of two component fields. One of the 
fields arises due to the inhomogeneous density of the 
condensate and leads to the precession of an off center 
vortex around the trap center [23] [25]- In addition to 
this, each vortex has a velocity field which varies inversely 
with the distance from its center, which is experienced by 
the other vortex of vortex dipole. In the present work, 
we move the obstacle along a:-axis and generate vortex 
dipoles located symmetrically about ic-axis. If (x, y) and 
(x, —y) are the locations of the positively and negatively 
charged vortices of the vortex dipole, respectively, then 
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FIG. 1. (Color online) Stationary state |^»| of binary conden- 
sate with obstacle potential at (a) — 6.0a osc , (b)— 5.9a osc and 
(c) — 5.0a OS c. 



the velocity field of the positively charged vortex is 

v(x,y) = u px k x r + —i, 

where w pr is the rotational frequency of a vortex with 
charge +1 in the condensate. A similar equation de- 
scribes the velocity field of the negatively charged vortex. 



OBSTACLE MODIFIED DENSITY 

To examine the density perturbations from the ob- 
stacle beam, let R[ n be the radius of the inner species 
or the interface boundary. And, let i? ut be the radial 
extent of the outer species. In the absence of the ob- 
stacle beam, the chemical potential of first and second 
species in scaled units are fii — i?? n / 4 + u 11 /(7ri?? n ) and 
R 2 ut /2, respectively. The obstacle beam initially (t = 0) 
located at (— i? ou t,0) traverses towards the center with 
velocity w Q bs and the intensity is ramped down at the 
rate —dVo/dt = r\. The location of the beam at a later 
time is xo(t) = — i? out + Vo^t, and intensity of the beam 
is Vo(t) — Vo(0) — rjt, where Vo(0) is the initial inten- 
sity of the obstacle beam. At the starting point, the 
total potential V(R ou t, 0, 0) > R 2 nt /2 and the density 
of the outer species l^l 2 is zero around the center of 
the obstacle beam. However, as it traverses the con- 
densates with decreasing intensity, at some later time 
if, V(x (tf),0,0) < R 2 ou j2. Density |V> 2 | 2 is then finite 
within the obstacle. For compact notations, hereafter 
we drop the explicit notation of time dependence while 
writing xo(t) and Vo(t). 

A critical requirement to form coreless vortices is com- 
plete immersion of the obstacle beam within n\. Based 
on the previous discussions, as the beam approaches the 
origin, the last point of contact between the beam and 
interface at R m lies along x-axis. To determine the condi- 
tion when complete immersion occurs, consider the total 
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potential along x-axis around the obstacle potential 

r ,("o(i)) 2 



V(x,0,t) 



x 



Vo(t) 



1 



w- 



(x- Xo (t)r 



(3) 



where, the Gaussian beam potential is considered up to 
the second order term. The expression is appropriate in 
the neighborhood of the beam, and along ir-axis it has 
one local minima (x m i n ) and maxima each. There is also 
a global minima, however, it is not the correct solution as 
it lies in the domain where x > wj \/2 and hence outside 
the domain of validity of Eq. ^ . Correct global minima 
is located at x ~ and is associated with the harmonic 
potential. The obstacle is considered well immersed when 
x nlin is located at the interfacial radius i?j n , and let t; m 
be the time when it occurs. 

When the obstacle beam is well inside the inner species, 
within the obstacle beam n\ is zero but n 2 is nonzero. 
It then forms a second interface layer, which embeds a 
bubble of the n 2 within n\. Recollect, the first interface 
layer is located at R m and it is where n 2 encloses n\ . The 
second interface, unlike the one at R ln , is a deformed- 
ellipse and we label it as T. Around the interface, the 
two condensates mix with a penetration depth 
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and the density of the minority species decays exponen- 
tially. 

The transition from a single continuous interface to 
two separate boundaries at i?i n and T, when the obsta- 
cle crosses i?; n , is smooth in TF approximation and that 
is how we have defined t- lm . There are, however, strong 
perturbations when surface tension is considered, and the 
separation of the two interfaces occurs when the beam is 
deep inside n±. Prior to the separation, the interface is 
deformed to accommodate a long neck region where ri\ 
and n 2 are non zero. As the interface splits into two, 
there are large deformations from the equilibrium inter- 
face geometry, and surface tension generates a restoring 
force to bring it to equilibrium geometry. This creates 
density patterns with high curvature and initiates for- 
mation of the coreless vortex dipoles. 



Equilibrium TF n 2 within the obstacle potential at this 
instant of time is 
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This, however, is higher than the density distribution at 
t im , that is n 2r {x,y,t ini + At) > n 2r (x,y,t ini ) as the 
potential V is lower. This is on account of two factors: 
first, the amplitude of the obstacle potential decreases 
with time; and second, the harmonic oscillator potential 
is lower at Xo(t im + At). The number of atoms, however, 
does not change from the value at i; m unless there is a 
strong Josephson current. Density n 2 r is thus below the 
equilibrium value once the obstacle beam is well within 
m. This creates a stable bubble of n 2 assisted or trapped 
within the beam and is transported through the nx- 

Departure of n 2 from the equilibrium is not the only 
density evolution within the beam. There is a progressive 
change of nxr (density of first species within the obstacle 
beam) as the beam moves deeper into ri\. At time ti m , 
when the obstacle is completely immersed in ri\ the ef- 
fective potential, experienced by nx, V{x, y, t- lm )+n 2 rux 2 
is larger than fix- So, n\ is zero within the beam. How- 
ever, if the rate of ramping r\ is such that at a later time 
V(x,y, tin, + At) + n 2 rui 2 < fix, while the beam is still 
within ni, there is a finite nx within the beam. Since 
«12 > V aiia 22 f° r the condensate, in TF approximation 
the bulk values of nxr and n 2 r can not be simultaneously 
non-zero. At the same time, n 2 r is forbidden to migrate 
to the bulk n 2 due to the n\ generated potential barrier 
in the region between interfaces T and i?j n . To accom- 
modate both ni and n 2 within the beam, the shape of 
interface T is transformed to increase n 2 r- So that n 2 is 
zero in certain regions within the beam where the condi- 
tion, V(x,y,ti m + At) + n 2 Yii\ 2 < fix, is satisfied. This 
mechanism is responsible for obstacle assisted transport 
of n 2 across tlx- 



ENERGETIC STABILITY OF NORMAL VERSUS 
CORELESS VORTEX DIPOLES 



We use TF approximation to compare the energetic 
stabilities of normal and coreless vortex dipoles. 



OBSTACLE ASSISTED BUBBLE 



At a time At after the obstacle is immersed in m, the 
location and amplitude of the obstacle potential are 

x (t im + At) = -Ro nt + v b x (t im + At), 
V {t im + At) = Vb(0) - r) x (t im + At). 



Normal vortex dipole 

Assuming that the vortex affects the density of the 
condensate only within the core regions, we can adopt 
the following ansatz for binary condensate with a normal 
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where £ = l/^/2/ii is the coherence length of inner 
species. Using these ansatz the number of atoms are 
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In a similar way, we can evaluate the energy of the entire 
condensate. 



Coreless vortex dipole 



For coreless vortex dipole, we adopt the ansatz 
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X 2 + y 2 > Rf u 

[(x-v 1 ) 2 + (y±v 2 ) 2 ]^e 



FIG. 2. (Color online) \tpi\ and phase of binary condensates 
after the creation of coreless vortex dipole. (a) Coreless vortex 
dipole is fully formed but yet to dissociate from the obstacle, 
(b) center of obstacle beam is at origin and coreless vortex 
dipole is separated, (c) additional coreless vortex dipoles are 
generated when the obstacle reaches the interface, (d) phase of 
the binary condensate corresponding to (b), and (e) densities 
of the condensates parallel to x-axis and passing through the 
center of the coreless vortex. Blue arrow marks the center of 
the coreless vortex dipole. 



vortex dipole at (vx,±V2) 
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\x 2 +y 2 ^ R 2 n & 
\ {(x~v 1 ) 2 + (y±v 2 ) 2 ]>e, 
(R 2 n ^ (x 2 + y 2 ) < R 2 out || 
\[(x-v 1 ) 2 + (y±v 2 ) 2 ]^e 
(x 2 + y 2 )>R 2 ont (9) 
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Using these ansatz the modified expressions for N 2 is 
N 2 
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(2Rf n nl + 8/X1M2 - 8A? n MiM2 + 8m?m2 - 1 
[Mi - 4«2Mi) • (10) 
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As done earlier, we can also calculate the total energy E 
of the system. The important change in E is the inclu- 
sion of interface interaction energy £7j nt . It arises from 
the interface interactions at the cores of the vortex and 
antivortex. Based on Ref. 27J, 
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where P is the pressure on the circumference of the cores 
and 



(x 2 + y 2 ) > R 2 out 
(x 2 +y 2 )<R 2 n . 
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The vortex dipole contributes mainly through the kinetic 
energy of 0i , which may be approximated with the value 
of single species condensate given in Ref. 26 



E vd = In 
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In both the case, i.e. with normal and coreless vortex 
dipoles, the energy can be minimized with the constraint 
of the fixed number of atoms and R ln as a minimization 
parameter. For the parameters set a without obstacle po- 
tential, the coreless vortex dipole has lower energy than 
the normal vortex dipole and is shown in Fig. [3] for the 
vortex dipole located at (0, ±1). 
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FIG. 3. (Color online) The energy of the binary condensate 
with Vo — and rest of the parameters same as those in the 
parameters set a as a function of R[ n . The condensate has a 
vortex dipole located at (0, ±1). Black and blue curves are 
for coreless and normal vortex dipoles respectively. 

For Ni = N 2 — 10 6 , a n = 51ao, 022 = 99ao and rest of 
the parameters same as in parameter set a, condensate 
with the normal vortex dipole has lower energy than the 
one with coreless vortex dipole (see Fig. |4| . These results 
are in very good agreement with the numerical results. 

66.6 1 . . . 1 . 1 
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FIG. 4. (Color online) The energy of the binary condensate 
with TVi = N 2 = 10 s , an = 51a , a 22 = 99a , V = 0, and 
rest of the parameters same as those in parameters set a as 
a function of R ln . The condensate has vortex dipole located 
at (0, ±1). Black and blue curves are for coreless and normal 
vortex dipoles respectively. 



NUMERICAL RESULTS AND CONCLUSIONS 

To examine the formation of coreless vortex dipoles 
in finer detail, we resort to numerical solution of Eq. ([2| 
with a modified version of the split-step Crank-Nicholson 



code reported in Ref. [22]. Consider obstacle potential 
is initially located in the outer component and is moved 
across the interface, towards the origin. For this case, 
we consider 85 Rb- 87 Rb binary condensate with param- 
eters set a, however, with maximum value of obstacle 
laser potential Vb(0) = 125.0. The obstacle potential is 
initially located at x = — 15a osc - The obstacle moves 
with the speed of 180/im/s, progressively decreases in 
strength with rate constant n — 10.1 (in scaled units), 
and vanishes at x = 8a osc . The obstacle potential cre- 
ates a normal vortex dipole as it traverses ni. As the 
obstacle penetrates the interface, it carries the vortex 
dipole generated in the outer component in its region of 
influence. Further motion of the obstacle, in m, creates 
coreless vortices. 

The key factor which influences the generation of core- 
less vortex dipoles is the deformation at the aft region of 
the obstacle confined 712 . The deformation accompanied 
by large mixing is initiated when the interface is about 
to break up. This is evident even in the stationary state 
density distribution shown in Fig. [ijb). At break up, 
the interface repulsion and potential gradient are highest 
along the x-axis and lead to the formation of a dimple. 
Curvature is large around the deformed interface, and 
flow of n% past it generates vortex dipoles. However, as 
the vortex dipoles are generated within the penetration 
zone, the build up of ri\ around the vortex core drives 712 
from the penetration zone to the core of the vortex. For 
the parameters considered, first coreless vortex dipole is 
formed soon after the interface break up and is shown 
in Fig. |2ja) at t = 0.21s. The vortex dipole is almost 
detached when the obstacle reaches origin, shown in Fig. 
[2jb). From the phase, Fig. [2jd), it is evident that the 
phase singularity is associated with n\ and ni is non-zero 
at the core. This is seen in the plot of densities along a 
line parallel to x-axis and passing through the vortex 
core, Fig. [2j». In the fi gure, the blue arrow marks the 
location of the vortex core. Another important density 
modification is, although the obstacle potential is repul- 
sive, ni has a maxima at the center. This is due to the 
repulsion energy from n\. More coreless vortex dipoles 
are created when the obstacle crosses the center of har- 
monic potential, [2jc) . 

The other initial configuration is to place the obstacle 
potential within n±. If V b s is such u\ < V(r,to) < ui, 
then within the obstacle potential ni is nonzero. Initial 
density distribution is like in Figs, [ljc) andJSJa). Due 
to inertia, ni lags behind the beam when the obstacle 
suddenly starts to move. When the obstacle has shifted 
from its initial position xq to xq + Sxq, the points of 
intersection of the inner interfaces (assumed circular with 
radii i?r and centered around Xq and x + 8xq) are 

(*«Ifc) = (l(2x +Sx),±^AR r 2 -Sx 2 ) . (13) 
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FIG. 5. (Color online) The generation of triply charged 
vortex dipole in the 85 Rb- 87 Rb binary condensate with pa- 
rameters set a. The obstacle potential, initially located at 
x — — 5.0a OS c, is moved with a velocity of 220^im/s up to 
x — 5a osc . First, second, and third columns are the solutions 
at t = Os, t = 0.14s and t = 0.32s, respectively. 



Due to inertia n 2 still occupies the region ST defined as 

- x Q ) 2 + y 2 < i?r < (x - x - Sx) 2 + y 2 

The potential experienced by ri2 along the left interface 
of ST decreases as one moves from x-axis. This leads to 
redistribution of ri2 in ST region and creation of a pres- 
sure difference 8P — (n 2 uu — n?,U22)/2 ^ at the left 
bounding arc of ST. Along this arc, SP decreases from 
(xq — Rr,0) to (x c ,y c ) and tends to flatten it. Another 
equally important dynamical process is the redistribution 
of nx as V e ff = V(x,y) + n 2 {x,y)ui2 increases along the 
left interface from x-axis to (x c , y c ). Due to this, n\ start 
to penetrate the interface from the point where V e f£ ^ Mi • 
Thus the repulsive interaction at the interface and gra- 
dient of harmonic potential, combine to form a dimple. 
The dimple formation initiates the formation of (coreless) 
ghost vortices in the obstacle region [35] , which detach to 
form coreless vortex dipoles as is shown in Figs. [5jb)-(c). 

We have studied the motion of the Gaussian obstacle 
across a phase-separated binary condensate. With the 
possibility of tuning one of the scattering lengths using 
Feshbach resonances, these condensates can be used to 
experimentally realise the obstacle assisted transport of 
one species across another as well as coreless vortices. 
Using both TF approximation and exact numerical solu- 
tions of coupled GP equations, we have shown that core- 
less vortex dipoles can be energetically more preferable 
than normal vortex dipoles. 

We thank S. A. Silotri, B. K. Mani, and S. Chattopad- 
hyay for very useful discussions. The numerical computa- 
tions reported in the paper were done on the 3 TFLOPs 
cluster at PRL. The work of PM forms a part of Depart- 
ment of Science and Technology (DST), Government of 
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